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ABSTRACT 

The observed X-ray source temperature distributions in OB stellar winds, as determined from high 
energy resolution Chandra observations, show that the highest temperatures occur near the star, and 
then steadily decrease outward through the wind. To explain this unexpected behavior, we propose a 
shock model concept that utilizes a well-known magnetic propulsion mechanism; the surface ejection of 
"diamagnetic plasmoids " into a diverging external magnetic field. This produces rapidly accelerating 
self-contained structures that plow through an ambient wind and form bow shocks that generate a 
range in X-ray temperatures determined by the plasmoid-wind relative velocities. The model free 
parameters are the plasmoid initial Alfven speed, the initial plasma-/? of the external medium, and 
the divergence rate of the external field. These are determined by fitting the predicted bow shock 
temperatures with the observed OB supergiant X-ray temperature distribution. We find that the 
initial external plasma-/? has a range between and 2, and the assumed radially-decreasing external 
magnetic field strength that scales as r~^ has a value of S lying between 2 and 3. Most importantly, 
the initial plasmoid Alfven speed is found to be well-constrained at a value of 0.6 Voo, which appears to 
represent an upper limit for all normal OB stars. This intriguing new limit on OB magnetic properties, 
as derived from Chandra observations, emphasizes the need for further studies of magnetic propulsion 
mechanisms in these stars. 

Subject headings: shock waves - stars: early-type - stars: magnetic fields - stars: winds, outflows — 
X-rays: stars - X-rays: general 



1. INTRODUCTION 

Waldron & Cassinelh (2007, 2008; hereafter, WC07 
and WC08) presented a detailed study of the Chandra 
HETGS data from 17 early-type stars (hereafter, OB 
stars). One of the more surprising results is the evi- 
dence for a radially dependent X-ray temperature (Tx) 
distribution that steadily decreases outward through the 
stellar wind from a high value near the star. We suggest 
that this distribution may be related to stellar magnetic 
field effects, and present a semi-empirical model for X-ray 
production in OB stars that incorporates a well known 
magnetic ejection process. 

The plausibility of magnetic ejection is supported by 
the increasing number of magnetic field detections on 
OB stars. There are now direct measurements of mag- 
netic fields for three 0-stars (Donati et al. 2002, 2006a; 
Bouret et al. 2008) and five early B-stars (Donati et al. 
2001, 2006b; Neiner et al. 2003; Petit et al. 2008). The 
recent observation by Bouret et al. (2008) provides the 
first detection of a magnetic field on an O supergiant, 
C Ori A (09.7 lb), with a field strength of 60 - 100 G. It 
also has a complex magnetic topology apparently simi- 
lar to the magnetic field structure of r Sco (BO V) that 
is interspersed with "hole regions" that may resemble 
the magnetospheric structure of the Sun (Donati et al. 
2006b). This C Ori A detection is of particular interest 
to us since magnetic field strengths similar to the ob- 
served value were proposed to explain its high energy 
X-ray emission lines in Einstein SSS data (Cassinelli & 
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Swank 1983) and Chandra HETGS data (Waldron & 
Cassinelh 2001). 

Theoretical studies of magnetic fields in OB stars 
have primarily focused on studying origins; dynamo pro- 
cesses (Charbonneau & MacGregor 2001; MacGregor 
& Cassinelli 2003; MuUan & MacDonald 2005), and 
remnant fossil fields (Ferrario & Wickramasinghe 2005, 
2006). In addition, Maeder et al. (2008) has shown that 
rapid rotation can extend the thickness of the thin outer 
convection zones in OB stars (Maeder 1980). At this 
early stage of studying OB surface magnetic field struc- 
tures we can only surmise that these field structures 
may likely be analogous to solar emerging magnetic field 
regions (Zwaan 1985). Unlike recent research in solar 
physics our intent is not to explore ways of producing 
coronal regions around these stars, but rather to explain 
the observationally derived X-ray temperatures of hot 
shock fragments that are dispersed throughout their stel- 
lar winds. 

2. THE Tx DISTRIBUTION IN OB STARS 

At the core of our semi-empirical model is the observed 
radial decreasing source Tx distribution found by WC07 
(and updated in WC08) which is most noticeable for OB 
supergiants (see Fig. [T]) . This distribution is unexpected 
based on wind shock theory developed prior to the launch 
of Chandra (e.g., Feldmeier 1995). To analyze the sig- 
nificance of this Tx distribution we first obtain an em- 
pirical power-law fit (error-weighted) to the WC08 data 
(neglecting the S Ori A Ne IX data for reasons discussed 
by WC07) given by 

/ n \ 0.79±0.05 

Tx{r) - [18 ±2) (^] MK (1) 
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Fig. 1. — Comparison of the OB supergiant power-law fit given 
by Eq. {TJ (black line) with the observed Tx for all OB supergiants 
(SG), giants (G), and main sequence (MS) stars as determined by 
WC08. The two data points significantly outside this bound are 
associated with the peculiar MS star 9^ Ori C which has a very 
strong magnetic field of Ri 1000 Gauss (Donati et al. 2002). 

This power-law fit shows that very high Tx are present 
near the star and we think this expression will be use- 
ful for interpretations of the emission line results and 
for general modeling efforts. As shown in Figure [1] the 
lower luminosity class OB stars have more complex Tx 
distributions since these stars have winds that are more 
optically thin to X-rays, and there can be radially ex- 
tended emission regions associated with any Tx value. 
Nevertheless, Figure [T] shows that nearly all giant and 
main sequence stars have an upper bound on Tx, and 
this upper bound is consistent with the supergiant power- 
law fit (Eq. [T]). As discussed by WC07, clarifications of 
these observed Tx distributions should lead to a better 
understanding of the overall X-ray emission process in 
OB stars, and in § 3 we present a model that provides a 
viable explanation of this radially-dependent Tx upper 
bound. 

3. PLASMOID MODEL 

The primary goal of this paper is to find a process 
that can generate high X-ray temperatures at small stel- 
lar radii (< 1.5 i?*). In the wind-clump model of Howk 
et al. (2000) they realized a need for high Tx at small 
radii and chose to consider the possibility of in-fall. Here 
we consider an alternative method of producing large rel- 
ative velocities at small radii by using an early magnetic 
propulsion model suggested by Cargill & Pneuman (1984; 
hereafter CP84). They derived, in a straightforward way, 
the dynamics of isolated "diamagnetic plasmoids" that 
are launched into a diverging external magnetic field. 
The plasmoids are produced by magnetic reconnection 
events and are rapidly accelerated owing to the interac- 
tion of the plasmoid magnetic moment with the exter- 
nal field (Schluter 1957; Parker 1957; Pneuman 1983). 
We picture these plasmoids as being produced just above 
the photosphere, and due to the rapid acceleration, plow 
through the overlying stellar wind that has drawn out 
the external magnetic field. The plasmoid-wind interac- 
tions produce "bow shock" structures around the plas- 
moids that are likely to be similar to those described by 



Cassinelh et al. (2008; hereafter CAS08), and the Tx will 
depend on the relative velocities between these plasmoids 
and the wind. An advantage of bow shock produced X- 
ray emission is that the emission measure is determined 
by the influx of the ambient wind material. 

We neglect the effects of the line radiation force on 
these plasmoids because Abbott (1982) found that the 
line driving is dominated by those spectral lines that have 
energies slightly larger than the energy associated with 
the peak of the radiation field. Hence, for a parcel of 
gas that is in a higher stage of ionization than expected 
from the radiation temperature, the line force will be 
greatly reduced. The gas temperature of the plasmoid 
is likely to be higher than the surroundings because the 
magnetic reconnection process should generate sufficient 
energy to increase the temperature ionization state of the 
plasmoid (e.g., a gas temperature > 0.2 MK should be 
sufficient to shift the dominant ions beyond the radiation 
field peak of OB stars). The plasmoid is expected to 
remain hot since it will heat up as it accelerates as shown 
by CP84. In addition, the X-rays from the bow shock 
formed on the outer plasmoid surface will also contribute 
to maintaining the higher level of ionization by the Auger 
effect as discussed CAS08. 

The temperature of the X-ray source behind a shock 
depends on the relative velocity of the incident mate- 
rial and the shock front. Consequently, to produce the 
high Tx at low radii (< 1.5 i?*) where the ambient 
wind velocity is low, these plasmoids must experience 
very rapid accelerations to generate the required large 
relative velocities. As shown by CP84, the largest ac- 
celerations can be achieved only for the case when the 
initial plasma- /3 of the plasmoid, /3poj is << 1, where 

[3po = 87TPp{ro) / Bp{ro)^ is the ratio of the initial plas- 
moid gas pressure to the initial magnetic pressure at the 
radius (r^) where the plasmoid is created. Using the mo- 
mentum equation and the basic assumptions given by 
CP84, along with this condition on /3po, one can obtain 
an analytic expression for the plasmoid velocity (assum- 
ing the plasmoid velocity = at Tq), 

1/2 

where Vapo = Bp{ro)/{^iTPp{ro)) is the initial plas- 
moid Alfven speed, Pp{ro) = initial plasmoid mass den- 
sity, Ves is the stellar effective escape speed, and S is 
the exponent of the assumed radially-dependent exter- 
nal magnetic field, Be{r) — Be{ro) {ro/r) ■ 

Equation ([2]) is more general than the one given by 
CP84 due to the presence of 7(r) which arise from the ex- 
ternal gas pressure gradient contribution to the plasmoid 
acceleration. By allowing for a non-zero (3eo we derived 
equation ^ by starting with CP84 equation (2.16), and 
found that 7(r) is 



7(r) = 



1 



1 



1/4 



(3) 



where Pdr) = 8TrPe{r)/ Be{r) is the radially-dependent 
external plasma-/?, and (3eo = f^ei^o)- In the limit of 
Pe{i") (i.e., Peo and 7 1), equation ^ re- 
duces to equation (3.2) of CP84. In this initial study we 
assume that the external temperature is isothermal and 
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Fig. 2. — Parameter study of the plasmoid model Vp (top 
panel) and Tx (bottom panel) dependence on Vapo/Kx; = 
0.4, 0.5, 0.6, 0.7 for a fixed value of /3eo = 1. For each Vkpo/Voo, 
two curves are presented, for S = 2 (solid line) and 5 = 3 (dashed 
line). In top panel, the ambient wind velocity law is shown as a 
solid black fine. In bottom panel, the solid black line represents 
the power-law fit given by Eq. STl . 

equal to the stellar effective temperature and the over- 
lying wind is spherical symmetric which are common as- 
sumptions used in the development of basic wind models 
(e.g., Lamers & Cassinelli 1999). Hence, for Pe{r) > 0, 
7(r) will be < 1 and the plasmoid will experience an 
enhanced acceleration as discussed below (see Fig. 

Although bow shock structures generate a range of Tx 
around a plasmoid at any given radius (see CAS08), in 
this paper we are only concerned with the maximum Tx 
that occurs at the apex of the bow shock. From the 
Rankine-Hugoniot temperature relation, this maximum 
bow shock Tx(r) is 

Tx (r) = u( ^YM.\ MK (4) 

^ ' V 1000 J ^ ' 

where AV{r) — Vp{r) — Vw{r) (measured in kms~-^) 
is the maximum (apex) relative velocity between the 
plasmoid and the wind. For the ambient wind ve- 
locity (Vw) we adopt the commonly used expression, 
Vwir) = Foo(l - R*/r). For AV > 0, the plasmoid 
is moving faster than the ambient wind and the bow 
shock is formed on the outward-facing side of the plas- 
moid. Conversely, for AV < 0, the bow shock is formed 
on the star-ward side of the plasmoid since it is moving 
slower than the wind, a condition more likely at large 
wind radii. 

We present a series of models to illustrate the depen- 
dence of Tx on the three basic parameters {Vapo, /3eo, 
S). We use a T4o = 2158 kms'^ and a Ves = 782 
kms^^ which represent averages of the OB supergiants 
with measured Tx values. We set ro = 1.005 i?, which 



Fig. 3. — Same as Figure 2 but for the parameter study of 
Peo = 0, 0.5, 1, 4, using a fixed value of Vyipo/Voo = 0.6. The 
radius axis has been re-scaled to magnify the details at small radii 
where f3eo effects are most prominent. 

corresponds to an initial ambient wind velocity of 1/2 
times the thermal speed using a mean effective temper- 
ature of 36,000 K. The functional radial dependencies of 
Tx on Vapo/Voo, Peo, and S are shown in Figures [2] and 
[21 As is evident in Figure [21 the ratio Vapo/Voo is the 
dominant parameter determining the overall magnitude 
of the observed Tx distribution. Figure [31 shows that 
large changes in /3eo are required to affect the Tx distri- 
bution, but only below 2 An increase in /3eo produces 
a rise in the peak Tx and shifts it to smaller radii. The 
field divergence parameter, 5, primarily produces small 
vertical shifts in the Tx distribution. 

Non-linear least-squared fit models (A & B) are ob- 
tained by fitting the model Tx (Eq. [3]) to the ob- 
served OB supergiant Tx distribution using the super- 
giant Voo and Ves averages. The model fits are shown 
in Figure [H and best-fit parameters are listed in Table 
1. Model A uses all the data, and Model B neglects 
the very high Tx (« 26 MK) at « l.li?* that is as- 
sociated with Cyg 0B2 No. 8a. As evident in our pa- 
rameter study, both models predict a well-constrained 
value for Vapo/Vqo, emphasizing the importance of this 
parameter in explaining the Tx distribution. As for the 
other two parameters, the constraints are less stringent 
as shown in Table 1. The high f3eo is necessary to fit the 
Cyg 0B2 No. 8a low radius, high Tx value, as demon- 
strated in our /3eo parameter study shown in Figure [3l 
As shown in Figure [H the two model fits (A and B) are 
essentially identical above 2 i?, , whereas, below 2 i?,, 
the models are more sensitivity to the values of /3eo and 
S. 

From the derived value of Vapo = 0.60 Voc we can check 
the assumption that Ppo << 1, which is fundamental to 
the validity of equation ([2]), and we also derive a plasmoid 
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TABLE 1 

Plasmoid Model Best-Fit Parameters 



Model 




/3eo 


S 


xVdof 


Model A 


0.61 ± 0.01 


1.3 ±0.9 


2.2 ±0.3 


17/13 


Model B 


0.60 ±0.01 


0.28 ± 0.31 


2.6 ±0.2 


15/12 
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Fig. 4. — Comparison of the two plasmoid ejection model fits to 
the observed Tx distribution of OB supergiants; Model A (solid 
blue line) and Model B (dashed blue line). See Table 1 and discus- 
sion in §3. The solid black line represents the power-law fit given 
by Eq. Q. 

initial density. By writing /3po — '^{ap/VApo)'^ — 0.02 Tg, 
where Up is the plasmoid initial thermal speed (Tg has 
units of MK), a value of (3po < 0.1 requires a < 5. To 
satisfy the requirement for an enhanced ionization state 
in the plasmoid as discussed earlier, all that is needed is 
a plasmoid temperature > 0.2 MK, well below this up- 
per limit. Assuming that the detected field of « 80 G 
from ^ Ori A is typical for the initial external field in all 
OB supergiants, we estimate from Vapo an OB super- 
giant plasmoid initial density of 4 x 10^°(1 + (3eo) cm~^ 
where we have used the total pressure balance equation, 

Bp{ro)^ — Be{ro)^{i + Peo)- Since this density is much 
lower than that expected for the surrounding external 



gas close to the star, these plasmoids should not be vi- 
sualized as clumps, but more appropriately, as isolated 
magnetic rarefactions. 

4. DISCUSSION 

We have presented a novel idea that can explain 
the near-surface high X-ray temperatures [i.e., NSHIP 
( "near-star-high-ion- problem" ) introduced by WC07], 
and the overall radially-decreasing Tx distribution ob- 
served in OB supergiants as deduced by WC07. By 
using a mechanism suggested by CP84, self-contained 
diamagnetic plasmoids produced near OB stellar sur- 
faces can be rapidly accelerated to speeds in excess of 
the ambient wind terminal velocity within 2 i?*. These 
high speeds can only occur in the CP84 theory when 
the ejected plasmoids are initially magnetically domi- 
nated (i.e., (}po « 1). The acceleration is most strongly 
dependent on the plasmoid initial Alfven speed (Vapo)- 
Since our model fits demonstrate that a specific limit of 
Vapo < O.6V00 provides a physical explanation of the ob- 
served radially-dependent Tx upper bound for all normal 
OB stars, it seems that we have found a new upper limit 
to the initial plasmoid Alfven speed for these stars. In 
addition, our derived Peo and S parameters provide lim- 
its on the initial external plasma conditions. Our results 
demonstrate that Chandra observations have the poten- 
tial to provide interesting information on the plasmoids 
and external plasma properties, and provide magnetic 
field constraints which are particularly relevant for those 
OB stars for which Zeeman measurements are not yet 
possible. Given the successful aspects and predictions of 
this plasmoid ejection model, further observations and 
detailed dynamical studies of magnetic propulsion mech- 
anisms are clearly warranted. 
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